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Infrared and visible (IR-Vis) light beams carrying orbital angular momentum (OAM) 
are currently thoroughly studied for their extremely broad applicative prospects, among 
which quantum information, micro-machining and diagnostic tools. Here, we largely 
extend these prospects, presenting a comprehensive study for the synthesis and full 
characterization of optical vortices carrying OAM in the extreme ultraviolet (XUV) 
domain. We establish the upconversion rules of a femtosecond IR-Vis helically phased 
beam into its high order harmonics, showing that each harmonic order carries the total 
number of OAM units absorbed in the process. To demonstrate a typical use of these 
new XUV light beams, we show our ability to generate, through photoionization, and 
control attosecond electronic wavepackets carrying OAM. These unique attosecond light 
and electron springs are promising secondary sources for applications and fundamental 
tests. 


I. INTRODUCTION 

Like any massive particle may carry two types of an¬ 
gular momenta (AM), namely spin and orbital angular 
momenta (SAM and OAM respectively), the massless 
photon also possesses two types of angular momentum 
[1]. It was early recognized that SAM is associated to 
the circular polarization of light beams but, only lately 
was the OAM of light beams associated with tilted wave 
front beams and made largely accessible experimentally 
in the infrared and visible (IR-Vis) spectral regions [2]. 
Helically phased light beams with a phase singularity as¬ 
sociated to zero intensity along the axis are the most 
common implementations of macroscopic beams carry¬ 
ing OAM [3]. At the photon level, the twisted wavefront 
of the beam translates into a quantized momentum, with 
its component along the beam axis taking discrete ih val¬ 
ues, where £ is any positive or negative integer number 
[3]. The unique properties of twisted light beams which 
are showing up both at the macroscopic and photon lev¬ 
els, lead to the emergence of countless applications, rang¬ 
ing in the Vis-IR domain from quantum information [4] 
to microscopy [5], nanoparticles manipulation [3] or fine 
structuring of materials using pulsed lasers [6]. In the 
XUV domain, twisted light beams hold the potential for 
specific fundamental and applied prospects, for instance 
new kind of dichroisms [7, 8] or the visualisation of dis¬ 
location strain fields in bulk material [9]. These new 
perspectives have already motivated the development of 
XUV sources carrying OAM on both quasi-continuous 
synchrotron installations [10], and on femtosecond free 
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electron lasers [11, 12]. One of the main outcomes re¬ 
ported here is the observation of such beams both in the 
attosecond regime (1 as = 10“^^ s) and in the extreme 
ultraviolet (XUV) domain. 

As a table top, ultra-short and largely tunable al¬ 
ternative to large scale instruments, High-Harmonic 
Generation (HHG) based XUV sources show unrivalled 
stability specifications, especially useful for pump-probe 
experiments targeting attosecond time resolution [13]. 
They are based on the upconversion of a high intensity 
femtosecond Vis-IR laser beam into the XUV range 
through nonlinear interaction with a gas. Very recently, 
two groups reported HHG with helically phased light 
beams [14, 15]. However, the conclusions of these 
two papers as for the harmonics’ OAM are in strong 
disagreement. On the one hand, analysing a couple 
of low HHG orders, Gariepy et al. [15] confirmed the 
theory of Hernandez-Garcia el al. [16] and the general 
“multiplicative” rule for OAM transfer in nonlinear 
processes [3]: for the two harmonic orders they could 
characterize, they observed that the OAM of harmonic 
order q is the quantum number £q = g^i, where £i 
is the OAM of the driving field. On the other hand, 
Ziirch et al. [14] observed on a single harmonic that 
£q = £i^ arguing that parametric instabilities prevent 
higher order vortices from propagating towards the 
detector. XUV light pulses of attosecond duration 
carrying OAM, based on a broad comb of phase-locked 
helically-phased harmonics, thus remain to be observed 
experimentally. Such attosecond light springs [17], lying 
in the XUV spectral range, would be a unique source 
to tailor attosecond electron beams carrying OAM 
through photoionization. Twisted electron beams are 
actively studied for future applications in fields covering 
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spectroscopy of diluted and condensed-matter [18], 
microscopy and particle physics [19]. Currently, they 
are primarily generated through tailoring non-singular 
electron beams in the quasi-static regime [20, 21]. 

In the present work, we first report on the validity of 
the multiplicative 0AM transfer rule for HHG with in¬ 
frared helically phased beams over an extremely broad 
spectrum and in several conditions of practical interest. 
Thus, we establish the opportunity to synthesise XUV 
“light springs” with attosecond structures. We addition¬ 
ally use these light springs to generate attosecond elec¬ 
tron beams with vortices, which were fully characterized 
in time. Finally, we propose some promising applications 
for these beams such as tests of the photoionization se¬ 
lection rules as recently theoretically discussed [22, 23]. 


rules the order dependence of the harmonics rings radii. 
In particular, the multiplicative law £q = q x ii yields a 
constant ring diameter over the whole spectrum while 
(^1 7^ 0) would lead to a diameter increasing 
with q. It is worth noting that none of these predictions 
matches the behavior observed when driving HHG with 
regular Gaussian beams, which leads, for the central 
part of the XUV beam, to an increasing divergence with 
q. This regular behavior is mainly a consequence of the 
HHG process itself combined to the A-dependence of 
light propagation to the far field [24]. Here we predict 
that the helical phase would modify this regular behav¬ 
ior. Inspecting the divergence of harmonics should thus 
answer the controversy about the transfer law, without 
resorting to any diffractive element or interference 
scheme. 


II. MULTIPLICATIVE RULE FOR 0AM TRANSFER 
THROUGH HHG 

As raised up by Gariepy et al. [15], the main difficulty 
encountered when studying XUV twisted attosecond 
pulses lies in the need for 0AM characterization over 
a broad spectral range. In the Vis-IR range, the phase 
fronts of twisted light beams are usually analysed 
directly using diffractive optical elements or interference 
schemes to determine the amount of 0AM [3]. Ziirch 
et al. [14] and Gariepy et al. [15] transferred these 
approaches to the XUV domain, but soon faced the 
challenge of having to resolve fringes with extremely 
low spacing due to the short wavelength of the XUV. 
In practice, the 0AM of only one or two low HHG 
orders were analysed. However, these light beams show 
ring-shaped intensity profiles while propagating, which 
is another signature of their 0AM. Here, we exploit 
the parameters of these peculiar profiles for their 0AM 
measurements. 


To illustrate our approach, we focus on Laguerre- 
Gaussian (LG) modes, which are eigen solutions of the 
paraxial wave equation and form a natural set of helically 
phased light beams carrying all possible integer values of 
0AM [3]. We restrict ourselves to LG modes with a ra¬ 
dial index p = 0, thus showing a single intensity ring. 
With the simple hypothesis that HHG driven by such an 
IR LG beam with 0AM ii leads to a comb of XUV LG 
modes, with the q-th order carrying an 0AM, the radius 
of the XUV beam reads [24] 
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The above considerations rely on the hypothesis that 
a pure LG mode is generated through HHG. To provide 
more steady grounds to this assumption, we performed 
a numerical simulation of HHG in argon gas within the 
strong field approximation (SFA), including macroscopic 
propagation, in a similar way as in [25] (see [24]). The 
results of the calculations are reported in Figure I. As 
expected, HHG occurs along the high intensity ring 
of the generating IR beam. The full width at half 
maximum (FWHM) of the profile obtained by a cut 
along a radius of this ring (hereafter called the thickness 
of the ring) is a fraction of the thickness of the ring of 
the driving field, a consequence of the high non-linearity 
of HHG. As for the spatial phase, a typical example 
is given by the 15^^ harmonic order (HI5), showing a 
spiral running 15 times 27r along the ring. This implies 
that the helicity of HI5 is 15 times the helicity of the 
driving field and that HI5 photons carry 15 times the 
0AM of the fundamental frequency photons. This 
behavior is consistent with the multiplicative law of 
nonlinear optics, and was observed here for all computed 
harmonics, from HI I to H33. To mimic the experiment, 
we simulate the propagation of the XUV beam towards 
an observation plane located 80 cm downstream the gas 
cell. When reaching the detector the harmonic beams 
still display a ring shape, suggesting that generation 
at focus is dominated by the emission of a single LG 
mode. As predicted by Eq. (I) with £q = qx£i, the ring 
diameter is identical for all harmonic orders. It only 
depends on the 0AM value of the driving field, which 
affects the diameter of the ring. The latter increases by 
a factor of 1.4 \/(^ when doubling the helicity of the 

phase, consistently with the facts that the 0AM carried 
by the harmonics is doubled and that the divergence 
goes like a/|^i|, according to Eq. (I). 


where z is the distance along the optical axis from the 
focal point and wq is the equivalent of the waist for a 
Gaussian beam. Eq. (I) establishes that the value of £q 


Getting back to the simplest microscopic model of 
HHG, by which an electron tunnels out of an atom close 
to an extremum of the driving field, then acquires energy 
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Figure 1. Full calculation of high order harmonic profiles generated with a driving laser carrying OAM. 

(A) Intensity and phase transverse profiles of harmonic 15 at the exit of the medium. We consider a spatially and temporally 
incident Gaussian beam of 6.25 mm waist and 50 fs FWHM duration, focused by a 1 m focal length lens, in the middle of a 500 
//m wide (FWHM) Lorentzian argon jet. The equivalent Gaussian beam waist at focus wq is 40 //m, leading to an equivalent 
Rayleigh range zr of 6.5 mm. The maximum gas pressure is 10 mbar and the laser peak intensity at focus is 1.5 x lO^^W/cm^. 
The observation plane is located 80 cm downstream of the jet exit. (B) Diameter of the harmonic rings in the far field for 
— 1 (red circles) and — 1 (blue triangles). The data for = 0 (purple crosses) is the short trajectory divergence 
experimentally measured using a Gaussian beam and a 2 m focal length lens. Insets: corresponding intensity profiles for a 
selection of harmonics. 


by quivering in this field and finally recombines with its 
parent ion emitting its excess of energy as XUV radiation 
[13], it may be anticipated that HHG does not generate 
pure LG beams, as pointed out by Ziirch et al. (14). 
Indeed, in the full computation, extra faint rings appear 
on the outermost part of the beams. They are decreas- 
ingly divergent with harmonic order until they merge 
with the central one [24]. Such patterns were lately ana¬ 
lyzed for HHG with regular Gaussian beams =0) and 
attributed, depending on the generation conditions, to ei¬ 
ther diffractive effects. Maker fringes in the time domain 
[26] or interferences of the emission from different quan¬ 
tum paths during the excursion of the electron around 
the ionic core [27]. By numerically turning off these tra¬ 
jectories one at a time, we observed that the outer rings 
are mostly present in the contribution of the more diver¬ 
gent long one. Being in a loose focusing configuration, 
this allows us to attribute them mostly to diffraction ef¬ 
fects, rather than interferences between the two quantum 
paths, thus letting the main conclusion of this paragraph 
unchanged. 


III. ATTOSECOND XUV LIGHT SPRINGS 

The numerical results above were tested experi¬ 
mentally on the LUCA laser server in Saclay. The 
experimental setup is described in Figure 2 (A). The 
laser beam is spatially filtered to yield a very close- 
to-Gaussian input mode [28]. A spiral phase mask is 
inserted to impose a staircase phase profile on this 
incoming beam, which can either be 1, 2 or 3 times 27r 


per turn, using two successive masks. Focusing such 
a phase shaped beam leads to a distribution of light 
about the focal point very close to a LG mode carrying 
respectively = ±1,±2,±3 units of OAM, the sign 
being determined by the orientation of the phase mask 
[29]. This is illustrated by the intensity profile reported 
in Figure 2 (B), the characterization of the OAM using a 
diffraction technique being described in [24]. Harmonic 
spectra obtained using an IR driving field with =1, 
2, 3 are displayed in Figure 2 (C). For each harmonic 
order and all OAM values of ^i, we observe a clear 
ring shape pinched in the horizontal dimension by 
the dispersion of the grating. The divergence of the 
harmonics, seen in the vertical dimension, appears to be 
constant throughout the whole spectrum, in contrast to 
what is commonly observed with a Gaussian beam. We 
measured the average ring diameters to be 1.01 ± 0.02 
mm, 1.33 ± 0.04 mm and 1.61 ± 0.01 mm for = 1, 
2 and 3, respectively. This is in agreement within 5% 
with the dependency predicted by Eq. (1), which 

leads to ratios of \/2 1A and \/3 1.7 between the 

diameters. Note that among all transfer laws we tested; 
only the multiplicative law provides such a quantitative 
agreement with the observations [24]. This extends the 
conclusions of Gariepy et al. (15) to the whole harmonic 
spectrum^ opening the way to the synthesis of attosecond 
light pulses carrying extremely high values of OAM. 

Furthermore, it appears in Figure 2 (C) that the high¬ 
est observed harmonic order decreases with increasing £i. 
This drop in the cutoff energy is due to a lower IR peak 
intensity at focus, that can be attributed first to a lower 
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efficiency in the conversion from the IR Gaussian mode 
to the IR LG mode [29]. Second, the IR focal spot size 
increases with requiring a larger aperture of the IR 
beam to reach the nominal intensity for HHG. This is 
very demanding on the stability of the incoming phase 
profile of the IR field and, in practice, leads to imperfect 
intensity and phase profiles at focus. This translates into 
a lower yield of the XUV radiation and a degradation of 
the ring shape when increasing ii. However, the highest 
value of £q is obtained for £i=3 with harmonic 19, yield¬ 
ing £q=b7. Moreover, to get fairly high £q while keeping a 
smooth intensity profile, we observed that it is beneficial 
to rather use neon gas as a target with £i=l. The high 
ionization potential of this gas allows to couple a lot of 
energy into the target before its ionization is saturated 
while, due to the low ^i, it remains less sensitive to wave- 
front imperfections. As shown in [24], in these conditions 
we could generate harmonics up to the 41®^ order with a 
smooth intensity profile. Once again, we observed a con¬ 
stant divergence of the spectrum, confirming the general 
validity of our first measurements in argon. 


IV. ATTOSECOND ELECTRON SPRINGS 



Since the early days of attophysics it has been recog¬ 
nized that XUV attosecond pulses could tailor electron 
wave packets (EWP) through photoionization [13, 31]. 
Applications of such attosecond EWP were proposed, 
for instance, as a quantum stroboscope [32], or as tools 
to localize EWP in space and time around atoms and 
molecules [33]. In the present context, we propose to in¬ 
crease the tailoring knobs by transferring the 0AM car¬ 
ried by harmonics onto photoionized electrons, forming 
an electronic beam carrying an 0AM. The selection rules 
of photoionization when using a twisted beam carrying 
an 0AM have been theoretically investigated in differ¬ 
ent configurations in [22, 23] confirming that, for elec¬ 
trons out of the beam singularity, i.e. far from the focus, 
the electric dipole transition should remain dominant 
[24]. To test the validity of the electric dipole approx¬ 
imation, we performed an interferometric measurement 
of such electron beams, building on the well-established 
XUV-IR cross-correlation technique called RABBIT (re¬ 
construction of attosecond beatings by interferences of 
two-photon transitions) used to characterize attosecond 
pulse trains, see e.g. [30]. In RABBIT, a target gas is 
photoionized by the harmonic comb in the presence of 
a weak dressing IR field. The simultaneous absorption 
of an XUV harmonic photon (Hq) and one IR photon 
produces sidebands (SBq+i) in the electron spectrum at 
energies corresponding to even multiples of the driving 
laser frequency. The same sidebands are also populated 
by the absorption of an XUV photon from the follow¬ 
ing harmonic (Hq) and the stimulated emission of one IR 
photon. The interferences between these two quantum 


Figure 2. Experimental observation of HHG spectra 
carrying OAM. (A) Experimental setup. It has two op¬ 
eration modes, corresponding to the imaging of HHG light 
spectra and to the characterization of electronic wavepack- 
ets. The Ti: Sapphire laser beam is passed through a spiral 
phase plate (SPP) and used to generate harmonics. They can 
be either directly imaged using a photon spectrometer (Grat- 
ing-hMicro Channel Plate) or focused inside a magnetic bottle 
electron spectrometer time of flight (MBES). An extra weak 
IR beam may be superimposed with the XUV beam in the 
MBES with a delay controlled by a translation stage (TS), 
allowing us to perform RABBIT measurements. See [24] for 
more details. (B)Intensity profile of the laser beam at focus 
close to the HHG gas inlet. (C) Normalized intensity of har¬ 
monics 15th to 27th generated in argon and observed in the 
far field, using £i = 1 (top row), £i = 2 (middle row) and £i 
= 3 (bottom row). 


paths lead to oscillations of the sideband yield as a func¬ 
tion of the XUV/IR delay (tq) at twice the IR laser pe¬ 
riod. The phase of these oscillations is directly linked to 
the total relative phases of the two consecutive harmonics 
involved and to that of the two-photon dipole transition 
at play, dubbed as the "atomic phase". Most importantly, 
for the interferences to appear, the wavefronts of the IR 
and XUV need to be matched, such as a well-defined 
phase relation between the XUV and IR field across dif¬ 
ferent transverse parts of the gas jet is obtained. Using a 
twisted XUV photoionizing beam and assuming that the 
photoionized electron beam is also twisted (which would 
be a consequence of the validity of the electric dipole ap¬ 
proximation, see [24]), this condition would be met with 
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Figure 3. Attosecond electronic beams carrying OAM. (A) Two-color XUV+IR two-photon photoionization spectro¬ 
gram of argon for a driving field with The main lines correspond to odd harmonics while the weaker oscillating ones, 

showing a double periodicity of T= 1.33 fs and 2.7 fs are sidebands. The dynamical range of our spectrometer limits the number 
of measurable sidebands even if the harmonic cut-off lies higher. (B) Delay-averaged spectrum (red) in log scale and group 
delay of the electronic wave packets emitted (blue circles, the light blue strip about the curve represents the numerical analysis 
error bar at 3cr). The analysis of the sidebands is carried out with the theory of RABBIT [30], under the assumption of fully 
coherent light (34). (C) Spatio-temporal shape of the emitted electron beam using the intensity profiles from Fig. 2 and 
the group delays from panel (B). The 3D surface is a contour at 80% of the max intensity. The back panel is a projection 
of this intensity, (inset) Temporal cuts at three different azimuthal angles in the electron beam 0 = 0°(red), 40° (blue) and 
—60° (green). 


a dressing IR beam which is also twisted, but not with 
a regular Gaussian dressing beam. To be more specific, 
let us assume, as measured above, that the spatial phase 
of the q-th harmonic is (j)q{r^0) = where we use the 
(r, 6>)polar coordinates. The measured sideband inten¬ 
sity is an average of the contributions of each emitted 
electron in the interaction region. If we write the inten¬ 
sity coming from one point in this region, neglecting the 
contributions of the atomic phases we get [30]: 

SBg+i (cj, 0) = cos [2co’ro + ^q +2 ~ + {^q +2 — ^q — 2 ^ 1 ) 

where uj is the angular frequency of the driving laser 
and cpq the spectral phase of the qth order. To keep the 
intensity oscillating after integration over which is the 
operation mode of our detector, the 6>-dependent term 
must vanish, giving iq -^2 — ^q = 2^i. This condition is 
only met when taking iq = qii. 

We tested this prediction on the experimental setup 
based on a Mach-Zehnder interferometer schematized in 


Figure 2 (A). We could verify that even after reflection 
off a drilled mirror, the IR focus kept a donut shape and 
so did the harmonics, with again a constant divergence. 
Imaging the foci in the sensitive region of the electron 
spectrometer, we observed a thick donut profile for the 
dressing beam, ensuring a homogeneous dressing of the 
harmonics. A spectrum with sidebands lines was then 
observed (Figure 3 (A)). Scanning the XUV-IR delay, we 
observed 2uj oscillations of the sideband intensities. The 
observation of these oscillations supports i) that the XUV 
beam carries an OAM and has therefore a ring intensity 
pattern and ii) that the electric dipole approximation re¬ 
mains valid. As in standard RABBIT measurements, 
the phase of these oscillations is directly linked to the 
group delay (GD) of the EWP generated by the har¬ 
monics comb. This GD is plotted in Figure 3 (B). A 
linear dependence is obtained, corresponding to group 
delay between successive harmonics of Ate = 103 ± 9as. 
This value is almost equal to what has been reported 
with Gaussian beams {£i = 0) for such a driving peak 
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intensity [30]. 

The two sets of measurements reported above, i.e. the 
spatial intensity profile of the harmonics and their group 
delay finally provides us with all data needed to recon¬ 
struct the full spatio-temporal shape of the emitted at- 
tosecond electron beam. We hereafter focus on the case of 
photoelectron with linear momentum along the propaga¬ 
tion axis, copying the spring structure of the XUV. Here 
we suppose a flat cross section of argon, which is a good 
approximation in this energy range. As lately predicted, 
we get two intertwined spring like structures [17]. The 
helical shape is directly related to the ionizing XUV light 
helicity, while the presence of two such structures is rem¬ 
iniscent of the generation of odd harmonics only. A cut 
in the spatial domain displays a double lobe structure, 
while the temporal profile shows a pulse train with each 
pulse lasting about 200 as. An interesting characteristic 
of this structure is that it leads to an identical attosec- 
ond pulse train when observed at different azimuthal po¬ 
sitions, only delayed by half a period of the driving laser 
(1330 as) in 180 degrees, i.e. 7.4 as/degree. As shown 
in [16], the corkscrew structure still holds in the case of 
single attosecond pulses. We believe that this property 
makes such pulses powerful tool for transient absorption 
spectroscopy experiments, which require tunable pump- 
probe delays on the attosecond timescale. Here, one may 
spatially map the attosecond time-delay in a single shot, 
getting rid of any stability requirement. The dynamical 
range for the delay scan is here 1.33 fs, i.e. the infrared 
laser half optical cycle, which could easily be increased 
up to several femtoseconds by increasing the driving field 
wavelength. 


V. CONCLUSIONS AND OUTLOOK 

The experimental evidence presented here for the 
transfer of 0AM over an extremely broad spectral range 
down into the XUV region confirms recent theoretical 
predictions and ends ongoing controversies about this 
transfer law. In addition, the analysis method identi¬ 
fied here could easily be generalized to arbitrary large 
spectral bandwidth, opening the study of a variety of 
even less conventional experimental situations. Based on 
these grounds, of promising interest is the reported ob¬ 
servation of the attosecond structure of the generated 
XUV pulse train obtained through the measurement of 
its spectral phase. First, it confirms the conservation of 
the usual photoionization selection rules with helically 
phased XUV light beams. While we did not observe any 
specific behavior for the photoionization of a noble gas, 
whether or not the presence of 0AM might change the 
photoelectron emission times in more complex systems 
is unknown and deserves more theoretical investigations. 
Second, it opens the route to the manipulation of attosec¬ 
ond electron beams carrying OAM’s, which will use the 


large panel of attosecond tools developed in the past ten 
years, may it be through high harmonics spectroscopy or 
XUV-IR pump-probe experiments. In particular, specific 
dichroisms were predicted in the XUV spectral range and 
still remain to be observed. Multicolor HHG using such 
tailored beams could also provide all in one pump-probe 
schemes, taking advantage of the spatial encoding of the 
delay in the azimuthal phase with these beams. Finally, 
the precise characterization of these new light beams will 
pave the way to the study of the controversial coupling 
between SAM and 0AM in matter during photoioniza¬ 
tion. 
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Supplementary 

Materials 

MATERIALS AND METHODS 

The experimental setup is pictured in Figure 2 (A). 
The experiments were performed using the LUCA laser 
server at CEA Saclay which delivers 30 mJ, 50 fs, 800 
nm pulses at 20 Hz. The laser mode was converted 
using a 16-level spiral phase plate (SPP) manufactured 
by SILIOS Technologies. 

For RABBIT measurements, the laser is split in two 
uneven parts by a mirror with a 8 mm hole. The 
main (outer) part of the beam is focused by a Im focal 
length lens (L) into a gas jet provided by a piezoelectric 
driven valve (Attotech). A diaphragm (D) removes the 
remaining IR beam, while the harmonics are focused 
by a 0.5 m focal length toroidal mirror (TM) into the 
sensitive region of a 1 m long magnetic bottle electron 
spectrometer time of flight (MBES). A Si02 plate (SP) 
serves as further attenuation of the IR beam. An extra 
full weak IR beam (4 mm diameter, energy of 70 /iJ) may 
be superimposed with the XUV beam in the (MBES) 
with a delay controlled by a piezoelectric transducer 
(TS). 

For intensity measurements, the drilled mirrors are 
replaced by plain mirrors and L by a 2m focal length 
lens. The harmonics are collected downstream the 
MBES on a photon spectrometer made of a variable line 
spacing Hitachi grating (001-0437) and a micro channel 
plate (MCP) coupled to a phosphor screen. In order 
to reveal the 2D spatial profile of the harmonics while 
spectrally resolving them, the MCP are placed 8 cm be¬ 
fore the focal plane of the grating. The phosphor screen 
is imaged by a Easier Al02f CCD camera. The ob¬ 
servation distance from the source to the MCP is 115 cm. 


SI - DIRECT MEASUREMENT OF THE 0AM OF THE 
DRIVING LASER AT THE FOCAL SPOT 

The measurement of the 0AM in the visible domain 
usually relies on interferometric or diffractive schemes 
[3]. We here chose to check the behavior of our phase 
masks placing a triangular slit in the path of the beam. 
The experimental scheme used for HHG is just slightly 
modified: a triangular slit is placed before the focusing 
lens and the beam at focus is imaged on a CCD cam¬ 
era equipped with a microscope objective. As expected 
from the theory, when choosing a size of the variable slit 
comparable to the waist of the laser, we obtained a se¬ 
ries of spots along a triangular pattern, whose number is 
related to the 0AM carried by the beam (see Figure SI. 



Figure SI. Measurement of the OAM of the driving 
laser. The collimated and spatially filtered femtosecond 
laser beam propagates through one or two of the spiral phase 
plates (SPP). A variable equilateral triangular slit is placed 
about 10 cm downstream. The diffracted beam is then fo¬ 
cused by the same lens as for the HHG experiments. The 
focus is imaged on a Imagine Source CCD camera equipped 
with a xIO objective and a 160 mm long tube. The images 
recorded when placing the masks imposing a 27r, An and Ott 
azimuthal phases (corresponding respectively to^ = ^ = 2/i 

and i = Sh) are displayed as insets. In the latter case, two 
masks imposing 27r and An are twinned to yield a Ott phase 
per turn. 


This observation was only made when the incoming beam 
presented a very close to flat wavefront. For instance, 
constraining a mirror in its mount ruins this structure. 
This online diagnosis secured the quality of the incoming 
beam for HHG. 


S2 - 4D NUMERICAL SIMULATION LAYOUT 

In this paragraph we provide details about the simu¬ 
lation of HHG using a driving laser carrying an OAM. 
For the IR beam, we accounted for the full experimental 
setup described below and schematized in Figure 2. The 
structure of the LG mode, which does not possess the 
cylindrical symmetry of a Gaussian beam, required per¬ 
forming 4D (3D in space+time) simulations. To this end, 
the beam is propagated through the different optical el¬ 
ements, by means of the Huygens-Fresnel integral, up to 
the gas jet entrance. This is the starting point for HHG 
calculations. In brief the coupled propagation equations 
for the driving laser and harmonic fields are numerically 
solved on a three dimensional spatial grid, in the parax¬ 
ial and slowly varying envelope approximations, using a 
standard finite-difference method. Galculations are per¬ 
formed in a frame moving at the group velocity of the 
laser pulse, i.e. the speed of light in vacuum here. Both 
electron and atomic dispersions are taken into account. 
We first compute the driving field at a given time t in the 
pulse envelope and position z along the propagation axis. 
This field is then used for calculating the space- and time- 
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dependent ionization yields and dipole strengths. Ioniza¬ 
tion rates are modeled using Ammosov-Delone-Kra'inov 
(ADK) tunneling formula [34] while dipoles are computed 
in the SFA, following the model described in [35]. Once 
these two terms are obtained, the equation of propaga¬ 
tion for the harmonic field is solved. The calculation is 
repeated for each t and z. The spiral phase plate induc¬ 
ing the 0AM is modeled by (f{x, y) = ^itan“^^, where x 
and y are the coordinates in the plane perpendicular to 
the propagation axis. 


where = Ai/g is the wavelength of the q-th order. We 
suppose that LG IR driven HHG occurs at focus of a lens 
(z=0). Without giving any detail about the generation 
process, it can be assumed that the intensity profile of 
the XUV beam is maximal where the intensity of the IR 
field is maximal well. This is predicted in the Strong 
Field Approximation (SFA) [35] or directly by solving 
the Time Dependent Schrodinger Equation (TDSE) [36]. 
Equating the two maxima of the IR and the XUV yields 


S3 - ANALYTICAL EXPRESSION OF THE RING 
DIAMETER FOR AN 0AM BEAM 


w(Aq,0) 



w{Xi,0) 



Using cylindrical coordinates (r^O^z)^ where r is the 
distance from the optical axis, z the distance along this 
axis from the focal point and zr = t:wq/\i the Rayleigh 
range, wq the beam waist and Ai the wavelength of the 
radiation, the spatial profile of the electric field of a LG 
mode reads, up to normalization constants 


El (r, 6>, z) 


^iie 

w (A, z) 




Xi(z 


+ .2) 


^2 

w(X, 



(SI) 


where £ corresponds to the 0AM per photon in 
this mode, w{X,z) = -h {zj zr)^ is the equiva¬ 
lent of the waist for a Gaussian beam, and = 

(1^1 -h 1) {z/ Zr) is the equivalent of the Gouy phase. 
Let us consider HHG driven by such an infrared LG beam 
and assume that the q-th order is also a LG mode with 
an 0AM £q. The corresponding intensity writes 


Taking into account Xq = Ai/g, one gets: 


^max (^5 ^ d) 



/7rzre(Ai, 0)^ q£i 

\ ^ 


2 


which is the expression given in the text. We retrieve 
the fact that the 0AM of the beam rules its divergence. 
This is also obtained when using the generalized Hankel 
transform to calculate the propagation of the beam as 
detailed in the following section. 


Ie{r,e,z) 


w‘^ (A, z) \w (z) 


r ( 2r2 \ 


The maximum intensity along the radius is obtained at 
dll I dr = 0, i.e. 


4r 


w‘^ (z) J \w (z) 


2|q 


gV wZ(z) 




The solution to this equation is: rmax = (A, z) 

Eor each order, the maximum intensity is located at: 


'^max (Ag, z) — W (Ag, z) ^ 


= ^(Ag,0) 




1 + 


;(A„0)' 


S4 - PROPAGATION OF AN 0AM BEAM USING THE 
GENERALIZED HANKEL TRANSFORM 


In this section we describe how the propagation of 
XUV light carrying 0AM can be computed efficiently. 
With this treatment in hand, we can numerically com¬ 
pute harmonic profiles with all possible values of 0AM 
carried by each harmonic. We then find which solution 
for the “transfer law” matches our experimental data, 
which shows a constant ring diameter in the far-held 
over an arbitrary spectral range. 


The propagation of the harmonic radiation to the 
far held can be described using the standard Huygens- 
Eresnel integral. The electric held in a plane located at 
a distance z along the axis reads: 
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Ez, (x, y) = 




iXqZ 


oo oo 

f f A tt 2 I 2\ -i27r(xxQ+yyQ) 

/ / Eo{xo,yo) e ^3^'- “ ° -'e da;odyo, 


— OO —OO 


with (xo,^o) and {x^y) being the Cartesian coordinates 
respectively at focus and at the z position, Xq the wave¬ 
length of the considered harmonic, and Eq the harmonic 
electric field at focus. In the case of light carrying 0AM, 


this integral can be advantageously rewritten by consid¬ 
ering the symmetry of the problem. First we change 
to the cylindrical coordinates (i?, 0) for the far-held and 
(r, 0) at focus: 


Ez,\^ {R, 4’) 



iXqZ 


oo 271 

J J Eq e 

0 0 


— i2TTRr cos(0 —</)) 

zXq rdrdO 


In the case of HHG with a driving held carrying 0AM, 
symmetries may be invoked to write the XUV held at 
focus as Eq (r^O) = Uq (r) where Uq (r) is the r- 

dependent part of the complex held at focus. This yields: 






iXqZ 


oo 

I Uo (r) 


e 


271 


-i27TRr cos(6>-0) 

zXq rdrdO 


271 

Using the identity (x) = f cos a where 

is the n-th order Bessel function of the hrst kind [37], 
we can carry out the integration over 0 and obtain the 
intensity in the far held: 


^Z,Xq 0) 


Ott oo 

Y-fUo(r)e^ 

0 




/ 27TRr\ 



This integral now takes the well-known form of a Hankel 
transform. The (x) functions show aperiodic oscil¬ 
lations whose zeros location increases with ^q. In the 
following we track the origin of the constant divergence 
observed for beams carrying an 0AM down to the 
overlap of Uq (r) and ^ relevant Bessel 

functions are displayed in Fig. S2. 


For iq= 0, which is the usual case of HHG, the vari¬ 
ation of the extension of {R, (j)) with Xq is ruled 

both by the scaling of Jo ^ h Ihe phase 

terms in Uq (r). Uq (r) has a distribution of amplitude 
centered about r=0 in this case. Since Jo = 1 



Figure S2. Computation of the zeroth (A), hrst (B) and qth 
order (C-D) Bessel functions for arguments qx (A, B and D) 
and X (C). The x-scales have been adjusted to yield a few 
oscillations in each panel. The orders are q=13, 17, 21 and 
25. 


at R = 0 (Fig. S2 (A)), the main contribution to the 
integral thus come from this central part of Uq (r). 

However Jo varies little with the harmonic order 

(Fig. S2 (A)) and thus plays a minor role in the scaling 
with q. On the contrary, Uq (r) shows a large spherical 
phase term, whose curvature signihcantly varies with q. 
This is reminiscent of i) the driving laser intensity prohle 
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Figure S3. Computation of the Hankel transform in Eq. S2. The waist of the harmonic field at focus is taken to be 60 
mum. (A) Intensity of the 15-th harmonic during propagation along the optical axis. The field is propagated to the far-held 
using equation (S2) from z = 0 to z = 0.5 m. The role of 0AM in the held divergence is illustrated by taking either = 1 

(left) or -£i 5 = 15 (right). The intensity is normalized to one for every value of z and is plotted in logarithmic scale. 

(B) Relative diherence between the diameters of HI5 and H17 as a function of their respective 0AM. The far-held prohle of 
H15 and H17 are computed using (S 2 ) and their diameter are gathered for every value of ^15 and £±7 between 1 and 30. For 
every couple of (^ 15 ,^ 17 ) values, the diherence of their diameter is plotted in percentage of the diameter of H15. The dashed 
red curve is a linear ht of the combinations for which this diherence is minimal. It corresponds to £±7 = 1.13^15 — The 

same was done for all harmonics from q=13 to q=25, which yielded the general result -£g +2 = 


at focus, which is parabolic as a hrst approximation and 
ii) the dependence of the excursion of the electron in the 
continuum with the amplitude of the electric held, which 
has a hrst order term. Considering this phase, it is found 
that the beam diameter increases with q for the shortest 
electronic trajectories. 

Now, when considering £q ^ 0 and a ring intensity 
distribution for Uq (r), the result is completely diherent. 
Like in the case of = 0, within the thickness of the 
ring, the phase of Uq (r) is spherical. However, now 
the Bessel functions are zero at R = 0 and the main 
contributions to the integral come from outer parts of 
the intensity distribution. If £q = £i , the variation of 
the Bessel function with q remains signihcant (Fig. S 2 
(B)). One thus expects the ring diameter to vary 
with q. On the contrary, with £q = g x ^ 1 , and an 
argument proportional to q (Fig. S 2 (D)), the position 
and width of the first lobe is almost stationary with 
q. This is the origin of the constant ring diameter 
observed. Note that this is the result of an interplay 
between the scaling factor of the argument and the or¬ 
der of the Bessel function (compare Fig. S 2 (C) and (D)). 

To put this analysis on firmer grounds, we used 
Eq. refeq:stwo to propagate a given field at focus with 
a chosen 0AM value by carrying out a pth-order quasi¬ 
discrete Hankel transform [38]. We used the standard 
radial distribution for a LG mode as Uo(r). The result 
of the computation varying z is reported in Fig. S3 (A) 
for the 15-th harmonic when choosing either ^15 = 1 or 
^15 = 15. The effect of the 0AM carried by the beam 
is dramatic, suggesting that the measurement of the di¬ 


ameter of the rings is highly sensitive to this 0AM. We 
repeated this computation systematically for harmonic 
orders 13 to 25 and for = 1 to 30 for each harmonic, 
finally getting a 2D matrix of ring diameters indexed by 
harmonic numbers and 0AM values. The relative vari¬ 
ations of these values were compared to each other. An 
example of such a comparison for orders 15 and 17 is dis¬ 
played in Fig. S3 (B). Clearly, a line of minimal variation 
appears, which nicely fits the slope £17 • This 

was observed for all other couples of harmonic orders (q, 
q+ 2 ): it appears that the ring diameter stays the same 
only when ^^+2 = which yields for any harmonic 

order £q = q£i. This demonstrates that the multiplica¬ 
tive law for 0AM transfer through HHG is the only one 
matching our experimental observations. 


S5 - SPIRAL WAVEFRONTS AND ATTOSECOND LIGHT 
SPRINGS 

In this paragraph we intend to give physical insights 
on the multiplicative law of 0AM established in the 
main text and on the formation of light springs with a 
double spiral structure. 

The origin of the multiplicative law may be inferred 
analyzing the HHG process in the time domain. We con¬ 
sider the simplest microscopic model of HHG, by which 
an electron first tunnels out of an atom close to an ex¬ 
tremum of the driving field, second acquires energy by 
quivering in this field and finally recombines with its 
parent ion emitting its excess of energy as XUV ra- 
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diation [13]. The phase slip of the emission with re¬ 
spect to the driving laser extremum (denoted (/:^xuv) is 
firstly determined by the last two steps of the process, 
namely the time spent by the electron in the contin¬ 
uum and the phase of the recombination dipole. We 
denote T the laser period and temporarily consider an 
infinitely thin plane of radiating dipoles perpendicular 
to the main propagation axis of the laser beam. At dif¬ 
ferent azimuths due to the structure of the wavefront 
of the beam carrying an 0AM, the extrema of the field 
are shifted in time by ^i^T/27r (red stars in Figure S4 
(A)). The emission time of the XUV radiation (recombi¬ 
nation time of the electron) at azimuth 0 will thus write 
te{0) = {(fxuv^hO)T/27r^ showing a single spiral shape 
(green dots in Figure S4 (A)). In space, at a given time, 
the pitch of this spiral is Ai. To have an 0AM q x 
we should have q x ^isuch intertwined spirals within Ai. 
To retrieve this observation we need to consider i) that 
the medium is long compared to the wavelength and ii) 
that we are using a multi-cycle driving field (about 40 
extrema considering a temporal FWHM of 50 fs, as used 
in this work). If we consider emissions at both maxima 
and minima of all half cycles of the driving field when it 
propagates, the same sketch could be drawn T/q later. 
This delay is chosen so that at a given azimuth 0 the 
two emissions at time 0 and T/q later add up coherently. 
At T/q, the black curves representing the driving field, 
together with the red stars and green circles attached to 
them appear dephased by T/q while the blue curves re¬ 
main unchanged. This leads to a second dashed green 
line parallel to the first, offset hy Xi/q. As apparent on 
the right panel (Fig. S4 (B)), the final harmonic field 
thus shows a series of q intertwined spiral phase fronts, 
corresponding, as expected, to an 0AM q x ii. 

Finally, coherently superposing a series of successive 
odd harmonics having such a multiple spiral phase shape 
results in the destructive interference of most of these 
spirals, only two of them per cycle remaining (Fig S4 C 
and D). This is the origin of the intertwined double spiral 
shape displayed in Fig. 3 of the main text. 


S7 - EXPERIMENTAL SPECTRA USING NEON AS A 
GENERATING GAS 

High harmonics were generated in a neon gas target 
using ii = 1 for the generating beam. As seen in Fig. S5, 
we were able to generate harmonics up to 65 eV, corre¬ 
sponding to the 41-st order. The intensity exhibits a very 
nice ring profile, and we measure the average diameter to 
be 1.00 ± 0.05 mm, which is the same as what was mea¬ 
sured in argon. This confirms our former observations 
on a broader spectral range. The harmonics presenting a 
constant divergence throughout the spectrum, the mul¬ 
tiplicative law for 0AM transfer is verified. The highest 
harmonic therefore carries an 0AM of 41, while having 




z (units of X) 


Figure S4. (A-B) Plot of a series of progressively offset 
Acos z — iO curves against z at different azimuths 0 
going from 0 to 27r from top to bottom. Black curves: Ai, 
A=l, ^=1 and (j) = tt/S. Blue curves 5Ai, A=0.9, ^=5 and 
4> = 7r/3. On the left (resp right) panel, 7 (resp. 200) angles 
from ^ = 0 to ^ = 27r are represented. The red stars on the left 
panel are positioned at the maxima of the field, where tun¬ 
nel ionization may occur, and green circles are where recom¬ 
bination and XUV emission occurs. The right panel shows 
5 blue caustics, corresponding to the 5 intertwined spiral 
phase fronts. (C) Same as panel (A) but for 5 harmonics 
(3, 5Ai, 7Ai, 9Ai, llAi)) for each angle. (D) Color map of the 
intensity of the field resulting of the superposition of these 
five harmonics against the propagation axis (horizontal) and 
the azimuthal angle (vertical). 


a regular intensity profile. 



35 40 45 50 55 60 65 
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Figure S5. Normalized intensity of harmonics 23’^'^ to 41®^ 
generated in neon and observed in the far field, using i, — 1. 
The dashed white lines represent the average ring diameter. 


S8 - RABBIT WITH GAM’S 

Here we propose to analyze how the 0AM of a 
twisted photon beam transfers to electrons through 
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Figure S6. Geometry of the ionization of an atom by a twisted light beam. The coordinates (X,Y,Z) refer to the laboratory 
frame and (x,y,z) to the atomic frame centered on the atom nucleus. 


photoemission of atoms and investigate the validity of 
the dipolar selection rule. The formalism and notations 
below are closely related to the ones used in [22]. 

We first define the vectors which represent the electron 
and nucleus positions in the laboratory and the atomic 
frames, respectively. In the laboratory frame, the nu¬ 
cleus is located at R = (R, (a,/3) and the electron at 


r' = (r',^',0'), while in the atomic frame the nucleus is 
at (0, 0, 0) and the electron at f = (r, 0). The two sets 

of coordinates are related by r' = R + f. Rewriting Eq.(2) 
of [22] for p = 0 and using the vectors defined above, we 
obtain for the expression for the vector potential at the 
electron position r': 




1 

w (z') 



IGI 

exp 



Here, we have p' = r'sin^', z' = r'cos^', and the other 
quantities are the same as in Eq. SI. 


The photoelectron angular distribution for the absorp¬ 
tion of one photon is given by the differential cross sec¬ 
tion which, according to the Eermi Golden rule, requires 
computing such a matrix element 


J d^r A -P V’r(r), (S4) 


where (f), (f) are the initial and final state wave- 

functions, respectively. Two limiting cases appear: (i) 
Rsinof = 0; (ii) Rsina >> r la.u). In the first one, 
we recover the study of [22] with (p' = 0, which means 
that the angles appearing in Eq. S3 are the same as the 
ones in the integral of Eq. S4. Therefore, one obtains 
nondipolar selection rules (Eq. (6) of [22]) which depend 
on the value of the 0AM, reminiscent of the marked 
variations of A at the atomic scale, in the vicinity of the 
beam axis. In the other case, we can assume r' R, 
and retrieve the dipole approximation. There, A 
at the atomic scale becomes independent of the electron’s 
coordinates and just acts as a position-dependent ampli¬ 




p^z' 


^z’) Zr 


e') 


(S3) 


tude and phase factor, ^ ^Rj on the usual dipole ma¬ 
trix element. It leads to the standard AL = ±I rule, 
which is also used in different experiments involving sin¬ 
gle twisted photon transferring its 0AM to an ensemble 
of atoms [3]. This implies that, apart from a tight re¬ 
gion in space close to the laser propagation axis, where 
the electric field is very small, twisted photons act like 
"normal" (i.e. plane-wave) ones, except that they carry a 
R-dependent phase, which adds up in multiphoton pro¬ 
cesses like HHG or RABBIT. Thus, the q-th harmonic 
will have an 0AM of = q x and the R-dependent 
RABBIT oscillations of the sideband will be modulated 
by 


SBq+i (w, r) = cos [2 wTo + (y?q +2 - (/9q 

+ (V2-4-24)C(r)] . 

where is the 0AM of the dressing field and 
the usual harmonic phase. If the dressing beam is homo¬ 
geneous and shows almost fiat phase fronts {id = 0 for 
instance) the last phase term covers 0 — 27r and washes 
out the oscillations when averaging over the positions of 
the atom, as a temporal jitter would do. On the contrary. 
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dressing with an IR carrying an 0AM 


sures keeping the oscillatory behavior of the sidebands. 



